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Abstract Several exoplanets have been detected towards the Galactic bulge with the 
microlensing technique. We show that exoplanets in M31 may also be detected with 
the pixel-lensing method, if telescopes making high cadence observations of an ongoing 
microlensing event are used. Using a Monte Carlo approach we find that the mean mass 
for detectable planetary systems is about 2 Mj. However, even small mass exoplan- 
ets {Mp < 20 A/®) can cause significant deviations, which are observable with large 
telescopes. We reanalysed the POINT- AGAPE microlensing event PA-99-N2. First, we 
test the robustness of the binary lens conclusion for this light curve. Second, we show 
that for such long duration and bright microlensing events, the efficiency for finding 
planetary-like deviations is strongly enhanced with respect to that evaluated for all 
planetary detectable events. 
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1 Introduction 

Gravitational microlensing teciinique initially developed to search for MACHOs in the 
Galactic halo can be used to infer the presence of exoplanets around lens stars 
Indeed, the planet orbiting the primary lens star may induce significant perturbations 
to the standard (single lens) Paczyiiski like microlensing light curves [3l|3]. Until now, 
10 exoplanets has been detected towards the Galactic bulge (see http: / /exoplanet .eu") 
and the least massive planets have masses of about 3, 5 and 13 M®. The planet orbital 
separations are in the range 2-5 AU (about the Einstein ring radius). Microlensing 
technique complements the planet detections by other methods (Doppler radial velocity 
measurements and transits are among the most efficient) Q that are more sensitive to 
large planet masses (Jupiter-like planets) at small orbital distances. Microlensing gives 
also the opportunity to detect planets lying in M31 [5l[^. In this case, however, the 
source stars are not resolved by ground based telescopes, a situation referred to as 
"pixel-lensing" PTlf^EKS- Until now, however, only 25 microlensing events have been 
observed towards M31 [11] and in one case a deviation of the microlensing light curve 
from the standard Paczyriski shape has been observed and attributed to a secondary 
component orbiting the lens star [TS]. However, new observational campaigns towards 
M31 have been undertaken [13lll4lfT5] and hopefully new exoplanets might be detected 
in the future. 

The possibility to detect planets in pixel-lensing observations towards M31 has 
been already explored ^S alBul? . The analysis for planet detection, however, has been 
performed by using a fixed configuration of the underlying Paczyiiski light curve. In 
the present work, using a Monte Carlo (MC) approach [18] we explore chances to 
detect exoplanets in M31, by considering the multi-dimensional space of parameters 
for both lensing and planetary systems. By using the method of residuals, we select 
the light curves that show detectable deviations with respect to the Paczyhski shape. 
The advantage of the MC approach is that of allowing a detailed characterisation of 
the sample of microlensing events for which the planetary deviations are more likely 
to be detected. 

2 Generation of Planetary Microlensing Events 

Assuming that the lens is a planetary system lying in the M31 bulge or disk we generate 
a sample of pixel-lensing events towards M31. Microlensing parameters are selected 
according to the method outlined in [19II20I[2T| . To select planetary parameters (mass 
Mp and orbital period P) we adopt the relation [22] 

MMP,P) oc M,-P-^ {§) , (1) 

with a — 0.11 and /? = —0.27. This relation was obtained by investigating the dis- 
tribution of masses and orbital periods of 72 exoplanets (detected mainly by Doppler 
radial velocity method), and taking into account the selection effects caused by the 
limited velocity precision and duration of existing surveys 0. The upper limit of the 

^ At the moment about 400 exoplanets have been discovered, basically at relatively small 
distances from the Earth (see http://exoplanet.eu). 

^ A recent re-analysis 1231 of the mass-period relation employing a larger catalog of 175 
exoplanets detected by the Doppler radial velocity method is consistent with that of |22| . 
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planetary mass is set at 10 Afj which roughly corresponds to the lower mass limit for 
brown dwarfs. Moreover, following the numerical simulations for the planetary system 
formation [23j, we set the lower planetary mass limit at Mp =0.1 Mq. Once the 
masses of the binary components and the planet period have been selected, the binary 
separation dp is obtained by assuming a circular orbit. We also suppose that planetary 
perturbation time scales are much shorter than orbital periods of planets. All distances 
are normalised to the Einstein radius of the total mass A/ of the lens system 



where Dg (Dl) is the source (lens) distance. Moreover, it is assumed that all stars in 
M31 have planets. 

We simulate pixel-lensing light curves assuming a regular sampling (between 2 and 
24 hours), telescopes of different diameters in the range 2 — 8 m, typical observational 
conditions and an exposure time toxp of 30 and 60 min. Moreover, since in pixel- 
lensing towards M31 the bulk of the source stars are red giants (with large radii), we 
take into account the source finiteness by averaging the planetary magnification Ap{ti) 
(numerically evaluated by using the binary lens equation [251126] ) on the source size. A 
microlensing light curve is given by 

Sp{U) = hi + fo[{Ap{U))-l], (3) 

where /^^i is the background signal from M31 and the sky, /q is the unamplified source 
star flux and {Ap{ti)} the time varying magnification, that is averaged on the (pro- 
jected) source size E 

J^''deJ^'Ap{e,p-t,)i{p)pdp 

{Ap{t,)) = , (4) 

2^ Jo Hp)pd-P 

where p ~ 0^/9-^ is the normalised angular size of the source (Sg — Rg/Dg, i?g being 
the source radius, 9-^ = -Re/^L is the angular Einstein radius) and I{p) the intensity 
profile of the source including limb darkening (see [18] for more details). 

Planetary perturbations in the light curves occur when the source star trajectory 
in the lens-plane (the plane orthogonal to the line of sight to the M31 star source, 
passing at the lens position) crosses and/or passes near caustics. This is the caustic set 
of the source positions at which the magnification is infinite in the ideal case of a point 
source. Clearly for realistic sources of finite size the magnification gets still large, but 
finite. 

Light curves that show significant (> 3cr) flux variations with respect to the back- 
ground noise in at least three points identify a sample of "detectable" events. These 
curves are fitted with a standard Paczyhski form 

5°(i) = /w + /o[(^°(t.)>-l] , (5) 

where {A'^it.i)) is the usual magnification for a single lens and the averaging is made 
on the projected source size p. 

Next, to select light curves with detectable planetary features we consider two 
indicators: (i) the overall deviation (in units of error-bars) of the planetary light curves 
from the corresponding Paczyhski best fit, as measured by the sum of residuals Xr = 
X/i (*i)/-^tot (A'tot is the total number of points) and (u) the number of points 
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Fig. 1 Normalised (to unity) distributions of (top panel) and Rmax (bottom panel) 

for planetary selected events with Xr > 4, A^good > 3 ^i^d (e)max > 0.1. Here and in the 
following figures we select the following parameters: telescope diameter D=8 m, an exposure 
time ioxp = 30 min, a sampling time of the light curve of 2 hours. Solid lines correspond to I 
class events, dashed lines to II class events. 



■^goodi even not consecutive, which deviate significantly from the Paczyhski fit. Bj' a 
direct survey of many hgfit curves we fix the thresholds Xr th = 4 and N^^^^ q-^ = 3 and 
we select a sub-sample of light curves with detectable planetary features having Xr > 4 
and iVgood > 3. However, in this way, we find that some of the selected light curves 
show overall distortions (with respect to the Paczynski shape) which hardly could be 
attributed in real observations to a planet orbiting the primary lens star (we mean that 
it is not easy to reconstruct parameters of planetary systems from such observational 
data). This happens, in particular, for events with small planetary mass and/or large 
(projected) source radius, for which the finite source effects are maximised. A measure 
of the relevance of these effects is given by the relative (with respect to the Paczynski 
value) planetary magnification averaged on the source area 

//^d2^[|Ap(x,t,)-^0(x,t,)IM°(x,ti)]\ 
= [ J^. J ■ 

To be conservative, events with maximal finite source effects are eliminated from the 
following analysis by requiring, besides conditions (i) and (ii), that (in) (e)max > 0.1, 
where (e)max is the maximum value of (e(ii)) along the whole light curve. Conditions 
(i) and (ii) are particularly efficient to select light curves with a large number of points 
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Fig. 2 The maximum relative planetary magnification (e)max is given as a function of Mp. 



deviating from the Paczyhski best fit, tlie condition (in) ensures the presence on the 
hght curve of at least one clear planetary feature. 



3 Results 

In Fig. [1] (for D = 8 m), we give the distributions of ti^2 Rmax for the events with 
detectable planetary deviations (x? > 4, A'^good > 3 and (e)max > 0.1). Here, is the 
full- width half-maximum microlensing event duration ti^2 ~ *E ./(wq), where /(uq) is a 
function 9 of the dimensionless impact parameter uq and Rmax is the magnitude in the 
_R-band corresponding to the flux variation at the maximal Paczyriski magnification 
Af^max = foi^^max — !)■ In the Fig.[l]we discriminate two classes of events (indicated 
by I and II), according to the ratio p/uq > 1 (solid lines) or p/uq < 1 (dashed lines). 
The events of the I class with p/uq > 1 have short time durations {{ti/2) — 1-6 day) 
and larger flux variations ((J?max) — 20.6 mag). In these events planetary deviations 
are caused by the source trajectory crossing (in the lens plane) the central caustic 
region, close to the primary lens star. The events of the II class, with p/uq < 1, have 
longer durations {{ti/2) = 6.4 day) and smaller flux variations ((i?max) = 23.1 mag). 
Planetary perturbations in these cases are (mainly) caused by the intersection of the 
source trajectory with the planetary caustics and may also appear at times far from 
the maximum magnification time tg- 

The fraction of I class events is about 5%. However, it turns out that the prob- 
ability to have detectable planetary features in these events (that however are rare) 
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Fig. 3 Simulated I class event (points) and Paczyiiki fit (solid line). Some parameter values: 
p/no = 2.14, Rmax = 19.0 mag, ti/j = 0.6 day, Mp = 1.04 Mj, M = 0.55 Mq (M is the total 
mass of the lens system), dp = 3.6AU, i?E = 5.55AU, Xr = 46, (e)max = 0.2. 



is higher. This happens since the crossing of the central caustic is more probable in I 
class events with p/uq » 1. On the contrary, the generated events of the II class are 
more numerous, but have a smaller probability to show detectable planetary features. 

The existence of the two classes of planetary events in pixel- lensing towards M31 
is also evident in the Fig. [2] where we show (e)max (averaged on a large number of 
MC events of the same planetary mass Afp) as a function of Mp. We find that (e)max 
increases with increasing Mp, a result that is expected since the size of the caustic 
region is increasing with the planetary mass [2811291150] . We also verify that the finite 
size effects are more important for I class events. Indeed, for a given mass Mp, the 
events of the I class (with p/uq > 1) have smaller values of (e)max with respect to the 
corresponding II class events. 

An example of hght curve for a I class event is shown in the central panel of 
the Fig. O Here, we can see how the planetary deviations in the light curve (upper 
panel) obtained by using the point-like source approximation are washed out when the 
magnification is averaged on the source area (second panel). However, also in this case, 
planetary deviations, as measured by the residuals (bottom panel), still are present. 
The planetary event MOA-2007-BLG-400 shows a similar shape of the light curve [31] . 

An example of II class event is given in Fig. [4] where we see that also a small mass 
planet (Mp = 0.3 Afgj) can cause detectable planetary deviations in events for which 
the finite size effects are small (see also Fig. 8 in paper [3]; we remind that we have 
local light minima in the cases dp < 7?e)- In the considered event the geometry is such 
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Fig. 4 Simulated II class event (points) and Paczyhki fit (thin solid line). Some parameter 
values: pju^ = 0.04, Rmax = 24.0 mag, ti/2 = 18.7 day, Mp = 0.3 M^, M = 0.31 Mq (M is 
the total mass of the lens system), dp = 3.8AU, -Re = 3.3AU, Xr = 8.0, (e)max = 0.4. 



that the source trajectory is passing (in the lens plane) in the region between the two 
planetary caustics, where a deficit of magnification is present. Other examples of fight 
curves for I and II class events are given in [18| . 

The distributions of the planet mass Mp for the planetary selected events (x? > 4, 
-^good > 3 and (e)max > 0.1) is given in the Fig.[5](solid line). For comparison, the Mp 
distribution for the whole sample of detectable events (dashed line) is also given. From 
Fig. [5]it follows that larger planetary masses lead to higher probability for the detection 
of planetary features. This result refiects the fact that the detection probability is 
proportional to the caustic area (for small size sources), which increases with the planet- 
to-star mass ratio. From the same figure, it also follows that the exoplanet detection can 
occur with a non negligible probability for Mp > 0.06 Mj (Mp > 20 M®), although 
even Earth mass exoplanets might be in principle detectable. However, if we consider 
telescopes with smaller diameter {D < 4 m), the tail at low masses in Fig. [5] disappears 
and, practically no exoplanet detection occurs for Mp < 0.06 Mj. 

The probability of planet detection is maximised when the planet-to-star separation 
dp is inside the so called "lensing zone" , which is the range of star-to-planet separation 
0.6 < dp/i?E < 1-6 [271128] . The dp distribution for selected (solid line) and detectable 
(dashed line) events are shown in the upper panel of Fig. [6l The relevance of the lensing 
zone is clarified in the bottom panel of the same figure where the planet separation (in 
unit of the Einstein radius) d = dp/R^ is plotted. 
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Fig. 5 Normalised (to unity) distributions of the exoplanet mass Mp for the events with 
detectable planetary deviations (solid line) and for the generated events (dashed line). 



Our analysis shows that in pixel-lensing searches towards M31 the typical dura- 
tion of planetary perturbations is about 1.4 days. However, the number of significant 
planetary deviations on each light curve increases with increasing ratios p/uq. So, the 
overall time scale for planetary deviations can increase up to 3 - 4 days (for I class 
events). This means that a reasonable value of the time step for pixel-lensing obser- 
vations aiming to detect planets in M31 is a few (4-6) hours, almost irrespectively on 
the telescope diameter D. 

We also verify that the overall probability (we stress that here we are considering 
the full planetary mass range from 0.1 until 10 Mj, as well as the Rmax and t2/2 
ranges) to find pixel-lensing events with detectable planetary deviations is, however, 
very small: less than 2 % for Z) = 8 m and about 3%, if in the selection procedure we 
relax condition (in) (namely, (e)max > 0.1) and decreases rapidly for smaller telescopes. 



4 Analysis of the PA-99-N2 event case 

The POINT-AGAPE collaboration reported the detection of the microlensing event 
PA-99-N2 fSS^. This appears as a peculiar event first because of its extreme brightness 
[Rmax — 19 mag), long duration (ti/2 — 24 day) and location, some 22 arcmin 
away from the M31 center; second, as shown in a subsequent analysis [TJI, because the 

^ As discussed in Section 3, for the II class events the mean values of i?max and tif2 
about 23.1 mag and 6.4 day, respectively. 
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Fig. 6 Upper panel: normalised (to unity) distributions of the star-to-planet separation dp (in 
AU units) for the events with detectable planetary deviations (solid line) and for the generated 
events (dashed line). Bottom panel: distribution of dp//?E for events as before. 



anomaly with respect to the Paczyhki shape along the Hght curve can be attributed 
to a secondary component orbiting the lens star. In particular, An et al. [T2] have 
evaluated the a posteriori probability distribution for the lens mass which results to 
be extremely broad: for source and lens disc objects, they report a lens mass range 
0.02 - 3.6 Mq at 95% confidence level. Together with the small best fit binary lens 
mass ratio, q ~ 1.2 x 10~^, this puts the lens companion well in the substellar range. 
In paper [18] we have remarked that taking at face value the most likely value for 
the lens mass ~ 0.5 Mq for a disc lens, the lens companion would be a ~ 6.34 A/j 
object. This would make of the PA-99-N2 lens companion the first exoplanet discovered 
in M31. Furthermore, we had analysed the PA-99-N2 event within the framework of 
our simulation scheme, showing in particular that its (microlensing and planetary) 
parameters nicely fall in the expected range for II class events [18] • Here we further 
analyse this event. First, starting from the observational data0, we test the robustness 
of the binary-lens best fit solution. Second, we address the question of the efficiency 
for finding binary-like deviations for such bright and long duration events. 

In Table 1 of An et al. [12] a list of the best fit parameters adopting various 
models, including the best fit assuming a single lens, for different location of the lens 
system are given. In particular, the two best fit binary models, named CI and Wl, 



" Courtesy of the POINT- AGAPE collaboration. 
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Fig. 7 Normalised (to unity) distributions of the for binary (solid line) and single lens fit 
(dashed line). 



are indistinguishable in terms of their . An et al. [12] considered finite source 
corrections for the CI model (they named this model as FS) and found that these 
corrections decrease slightly the (~ 2%) with respect to CI model. Therefore, these 
authors claimed that CI and FS models are practically the same. In our analysis we 
consider the CI model. To verify the robustness of the binary-lens fit solution, the first 
test is to add a Gaussian noise to the best fit of the observed light curve to verify 
if a single lens model with noise can reproduce the observational data. We also let 
each parameter vary (by at most 20%) around the given central value. For each of the 
so obtained light curve we calculate the , which is plotted in Fig. [7] (dashed line). 
We have verified that the lowest value corresponds to the parameters of the best 
fit Paczynski model (last row) in Table 1 of An et al. [12]. Similarly, we have taken 
the best fit binary model CI and realised more than one thousand models by adding 
Gaussian noise and letting the parameters to vary by at most 20%. We plot again in 
Fig. [7]the distribution (solid line). Also here the lowest value is obtained for the 
CI model without noise. Moreover, from Fig. [7]one can see that the two (normalized) 
distributions of the are clearly separated, which implies that the best single lens 
fit is much worst than any of the binary lens models. From this we conclude that the 
binary fit is robust and that the observed light curve cannot be obtained by any single 
lens model with random noise. 



^ There exist several binary lens models that lie at a local minimum. Some of these 
degeneracy would have been removed if the ~ 20 day observation gap between JD' 72 and JD' 



91 had been regularly covered I12| . 
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Fig. 8 The binary light curves corresponding to the CI model [12J for g, r and i bands. 



A further point to be stressed is the following: at the end of the previous section 
we have mentioned that the probability to detect an exoplanet in pixel-lensing obser- 
vations towards M31 is rather small even with large telescopes. Therefore, the question 
which arises now is the following: what is the chance of finding a planetary feature in 
an event as PA-99-N2? The basic answer can be found looking at the characteristics 
of the underlying microlensing event. As it is also apparent by comparing with the 
distributions shown in Fig. [T] this is, at the same time, an extremely bright and long 
duration event. In fact, as we now show, this strongly enhances the probability for find- 
ing deviations to the single lens shape. To this purpose we perform a MC simulation 
where we fix the single lens parameters Rmax and ti^2 to those of PA-99-N2 and we let 
vary the binary ones (planet-to-star mass ratio, planetary distance in unit of -Re a-nd 
orientation of the source trajectory with respect to the binary axis) for an observational 
setup fixed to reproduce that of the POINT-AGAPE observations (mirror size, sam- 
pling and exposure time). In particular we find, for the selected planetary mass range, 
an increase of the average efficiency up to 6%, to be compared with less than 0.6% (see 
Table 3 in [T5]) for events without any constraints on -Rmax and ii/2- Moreover, as also 
implicit from the analysis of Fig. [5] the efficiency for finding binary-like deviations in 
the simulated light curves increases significantly with the lens-companion mass. In par- 
ticular, for the restricted range 1 — lOMj we find the average efficiency rises up to 27%. 
We have restricted our attention to the planetary-mass range for the lens companion 
(Mp < 10 Mj was our initial assumption), but we might expect the efficiency still rises 
once we enter the brown dwarf mass range (as in fact it is allowed by the analysis [12)'). 
Finally, the efficiency would be further enhanced for a lens-companion mass function 
more peaked towards higher values of the mass than that we have assumed. 
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5 Conclusions 

We have discussed the possibility to detect planets in M31 by using pixel-lensing ob- 
servations with telescopes of different sizes and observational strategies. Assuming that 
each lens star is hosting a planet, we carry out MC simulations and explore the multi- 
dimensional space of the parameters for both the lensing and planetary systems. Plan- 
etary induced perturbations in light curves are detected by comparing the simulated 
light curves with the corresponding Paczyhski shapes and searching for significant de- 
viations. The MC approach allows to characterise the sample of events for which the 
planet detections are more likely to be observed. Since in pixel-lensing towards M31 the 
bulk of the source stars are red giant (with large radii) , we take into account the finite 
source effects, which induce a smoothing of the planetary deviations and decrease thus 
the chance to detect planets. We estimate the typical duration of a single planetary 
feature to be of about one day. However, the number of significant planetary devia- 
tions and consequently the overall time scale of the perturbations increases (up to a few 
days), increasing the source size. Therefore for pixel-lensing searches towards M31 only 
few exposures per day could be enough to detect planetary features in light curves. 
Pixel-lensing technique favours the detection of large mass planets (Mp ~ 2 Mj), 
even if planets with mass less than 20 could be detected, although with small 
probability, by using large telescopes with a sufficient photometric stability. 

As a test case for exoplanetary searches in M31 we have reconsidered the POINT- 
AGAPE event PA-99-N2, which had already been probed to show an anomaly with 
respect to the Paczyriski shape compatible with a binary lens. As a first step we have 
revisited the issue of the single lens versus binary lens solution, finding that the lat- 
ter is indeed robust against the introduction of a gaussian noise along the observed 
data. According to the previous POINT-AGAPE analysis this binary system has a 
small mass ratio and this makes at least plausible that the lens companion is indeed 
an exoplanet. FHirthermore, the underlying microlensing event is extremely, and some- 
what unusually, long and bright. Therefore, as a second step, we have carried out a 
specific MC simulation that allowed us to show that for this kind of events the chance 
of finding exoplanetary deviations is indeed greatly enhanced and possible even for an 
observational set up as that of the POINT- AGAPE observations. As a caveat, we men- 
tion that the efficiency grows with the lens companion mass also beyond the exoplanet 
mass range in the brown dwarf one. 

Whatever the case for PA-99-N2 event, our analysis confirms that looking for exo- 
planets in M31 with pixel-lensing, at least in the Jupiter mass range, is already reach- 
able with present technology. Clearly, an efficient strategy of search, as towards the 
Galactic bulge, is mandatory: a wide field survey, to collect a large enough number of 
pixel-lensing candidates, endowed with an early warning system to trigger subsequent 
follow up observations, possibly with a network of telescopes around the world (for 
which also telescopes with small field of view could be usefully employed). In our opin- 
ion the reward of such a project would be substantial: going from the settling of the 
question of the MACHO fraction in M31 halo, to important information on the stellar 
mass function and the detection of cxoplancts, besides other information on the M31 
structure and content of variable stars. 
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